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Abstract 

We develop the plasmon-pole approximation for an interacting electron gas 
confined in a semiconductor quantum wire. We argue that the plasmon-pole 
approximation becomes a more accurate approach in quantum wire systems 
than in higher dimensional systems because of severe phase-space restric- 
tions on particle-hole excitations in one dimension. As examples, we use the 
plasmon-pole approximation to calculate the electron self-energy due to the 
Coulomb interaction and the hot-electron energy relaxation rate due to LO- 
phonon emission in GaAs quantum wires. We find that the plasmon-pole 
approximation works extremely well as compared with more complete many- 
body calculations. 
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I. INTRODUCTION 

Recently, there has been an increasing interestS in semiconductor quantum wire struc- 
tures, where the motion of electrons is essentially restricted to be one dimensional. Tech- 
nological progress has made it possible to fabricate! high quality quantum wires where only 
the lowest subband is populated by electrons, so that a truly one-dimensional interacting 
electron gas is realized. Much in the same way as quantum well structures have generated 
tremendous activities in pure and applied research on two-dimensional electron systems, 
quantum wire structures have created the potential for new device applicationJlli'i and the 
opportunity to carry out experimental study on one- dimensional Fermi systems, where many 
theoretical predictions! can be tested. Because of the low dimensionality, properties of a 
quantum wire are very sensitive to electron-electron interaction effectHi. Many experimen- 
tally relevant quantities need to be calculated by taking into account many-body interaction 
induced exchange-correlation effects. The standard perturbation theories, which have been 
developed for higher dimensional electron systems, have been appliedi'@ to quantum wire 
systems, and good agreement with experiments^'i are generally obtained. In this article, 
we discuss the application of another well known many-body approach, the plasmon-pole 
approximation,il0 to quantum wire systems. The motivation for this work is the observa- 
tion that the collective plasmon excitation plays a more prominent role in a one-dimensional 
electron system compared with its higher- dimensional counterparts because single-particle 
electron-hole excitation continuum is severely restricted in one dimension due to energy- 
momentum conservation. Thus, the plasmon-pole approximation, besides having the obvious 
benefit of great simplicity, may work well for one-dimensional quantum wires in calculating 
exchange- correlation effects. To illustrate this point, we calculate the electron self-energy 
correction due to electron-electron Coulomb interaction and the electron energy relaxation 
rate due to electron LO-phonon Frohlich interaction and compare the plasmon-pole approx- 
imation results with the corresponding full many-body calculations using the random-phase 
approximation (RPA). 
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The plasmon-pole approximation has been extensively employedBil3 in calculating the 
electron self-energies of three- and two-dimensional systems. The results obtained from these 
calculations are in good semiquantitative agreement with the results of more sophisticated 
treatments — namely the full RPA calculations, and with experimental results. A many- 
body interacting electron system has both collective plasmon excitations and single-particle 
electron-hole excitationslii. The plasmon-pole approximation simplifies the many-body ex- 
citation spectrum by ignoring the particle-hole excitations and assigning the whole spectral 
weight, which is dictated by the f-sum rule0, to an effective collective plasmon excitation, 
which is assumed to be a real pole of the response function. This is, in general, a crude 
approximation for the actual dynamical response of the electron system, except in the long 
wavelength limit where the plasmon excitation exhausts all the spectral weight in a uni- 
form systemlli by virtue of particle conservation. It is well known that in a one-dimensional 
system, the collective plasmon excitation plays a more prominent role because of phase 
space restriction on particle- hole excitations. In fact, the long wavelength RPA plasmon 
dispersion is exact in one- dimensional electron liquids up to second order in wavevector 
in contrast to higher-dimensional systemJl^. It is therefore worthwhile to explore the pos- 
sibility that the plasmon-pole approximation may actually work better in quantum wire 
systems than in higher dimensional systems. Our work is motivated by this purpose. We 
find that the one-dimensional phase-space restriction on the particle-hole excitations indeed 
increases the spectral weight of the plasmon excitation over a wide range of wavevectors in 
one-dimensional systems under conditions which are typical in GaAs-based quantum wire 
samples, and the plasmon-pole approximation indeed works extremely well in calculating 
a variety of quantities in GaAs-based quantum wire structures. We specifically apply the 
plasmon-pole approximation to two different problems: one involves the electron-electron 
Coulomb interaction and the other involves the polar electron LO-phonon Frohlich inter- 
action. In the first case, we calculate quasiparticle properties of the interacting electrons 
by taking into account the Coulomb interaction effects through the plasmon-pole approxi- 
mation. In the second case, we calculate the hot electron energy-loss rate via LO-phonon 
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emission. We show that the plasmon-pole approximation works well in both cases by giving 
results which are in good agreement with the corresponding RPA results, which are much 
more difficult computationally. 



system calculated within the RPA, where the particle-hole excitations are confined within 
the phase-space surrounded by the dotted-line ABODE, and the plasmon excitation is 
represented by the solid-line. In one dimension, the RPA plasmon dispersion has a simple 
analytical expression! (we set ^ = 1 throughout this paper) 



where E±{q) = q'^/2m±kFq/m with m as electron mass and kp the electron Fermi wavevec- 
tor, A{q) = exp[q / nVc{q)] with Vc{q) as the electron-electron Coulomb interaction potential^ 
of a quantum wire of finite lateral width. Unlike plasmon modes in higher dimensions, the 
RPA plasmon excitation in a quantum wire exists {i.e. is undamped) for all wavevectors 
< g < cx). 

The most characteristic feature of the one-dimensional spectrum is that particle-hole ex- 
citations are prohibited from a large portion of the low energy phase space, the region below 
the dotted-line BCD in Fig. |I[ This restriction, which arises from the momentum-energy 
conservation, increases the dominance of plasmon excitations in a quantum wire compared 
with higher-dimensional systems. The situation is totally different in higher dimensional 
systemJli, where particle-hole excitations are allowed in the whole phase space between the 
dotted-lines AB and DE in Fig. |l[ For a quantitative measure of the relative importance 
of the plasmon excitation, we evaluate its oscillator strength within the RPA 



where n is the average density of the electron gas, and the plasmon spectral weight is defined 
by (with the plasmon dispersion Ug given by Eq. (|l])) 



In Fig. 0, we show the elementary excitation spectrum of a one-dimensional electron 
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In Fig. 1^, we compare the oscillator strengths of plasmon excitation of a one-dimensional 
quantum wire with that of a two-dimensional quantum well (both obtained within the RPA). 
The plasmon oscillator strength in the quantum well drops quickly to zero at a critical 
wavevector, beyond which an undamped well-defined plasmon mode does not exist. The 
oscillator strength of the plasmon excitation in a quantum wire, on the other hand, extends 
well into the range of large wavevectors, decreasing slowly with increasing wavevector. Note 
that the q —* behavior of these curves, i.e. F{q) = 1 for g ^ 0, is just a manifesta- 
tion of the f-sum rule. The interesting point is that F{q) ~ 1 in one dimension even for 
q > kp. It is seen clearly that the plasmon dominance of the spectral weight is significantly 
increased in quantum wire systems. Since the plasmon-pole approximation assumes that 
the excitation spectrum consists of no particle-hole excitations, but solely of a collective 
mode which exists for all values of wavevectors and possesses unit oscillator strength, it is 
easy to understand why the plasmon-pole approximation may work well in a quantum wire 
system. The density- density response function of a quantum wire electron system in the 
plasmon-pole approximation is given as 

Xpp{q,uj) = ^^. (3) 

By construction, xpp(g, u;) in the above expression satisfies the f-sum rule and the static 
Kramers-Kronig relationlli. Eq. (|^) is our plasmon-pole approximation mo delHi, which we 
use to calculate quantum wire many-body electronic properties. 

In Sec. H and |ITTt we apply the plasmon-pole approximation to the calculations of 
electron self-energy due to Coulomb interaction and electron energy relaxation rate due to 
LO-phonon emission, respectively. A short summary in Sec. |IV| concludes our paper. 

II. SELF-ENERGY AND SPECTRAL FUNCTION 

The one-dimensional (ID) self-energy within the leading order GW approximationlll 
neglecting vertex correction at T = is given by 
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dqdoj' 
(2^ 



W{q,uj)GQ{k-q,uj -uj'), 



(4) 



where GQ{k,uj) is the Green's function for the noninteracting electron gas 

iki-kp) , e{kF-\k\) 



Go{k,uj) 
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(5) 



uj-^{k)+iO+ uj-^{k)-tO+' 

with C,{k) = k"^ /2m — /i, (/i=chemical potential) and W{q,uj) is the dynamically screened 
Coulomb interaction, which is given by 

Vc{q) 



W{q,u) 



(6) 



Here Vc{q) is the bare Coulomb interaction, which is logarithmically divergent in the ID 
wavevector space. Thus, we use the more realistic finite width quantum wire model whose 
fully approximated matrix element can be found in the literaturelll. e{q,uj) is the dielec- 
tric function, which describes the dynamical screening properties of the electron gas. The 
dynamically screened interaction W{q,uj) can be separated into an unscreened term which 
gives rise to the exchange part of the self-energy and another term which gives rise to the 
correlation part of the self-energy and involves coupling to density fluctuations 
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W{q,u) = VM) + VM) 



e{q,u) 



(7) 



The imaginary part of the second term is nonzero within the electron-hole continuum and 
along the plasmon dispersion line in the RPA. In the plasmon-pole approximation 
the second term is replaced by a coupling to the effective plasmon mode as described 
in Sec. | 



Im 



e{q,uj) 



--—6{uj-ujg), 

A UJq 



(8) 



where the strength Uq = ^Vc{q)q^ is determined by the requirement that Eq. (H) satisfies 
the /-sum rule, and Uq is the ID plasmon dispersion which is exactly known within RPA0 
(See Eq. (|I|)). Unlike in 2D and 3D, where the exact analytic RPA plasmon dispersion 
is unknown so that the static RPA dielectric function e(g, u; = 0) is used in obtaining the 



effective plasmon frequency uoq, we use the analytically known ID RPA plasmon dispersion 
given in Eq. (|T]). Note that any attempt to use the static RPA (similar to what is done 
in 2D and 3D PPA) in ID PPA is not only unnecessary (because the ID RPA plasmon 
dispersion is known analytically), but also incorrect because the ID static RPA dielectric 
function has logarithmic zero temperature singularities due to a divergence at g = 2kF- 
Using the Kramer-Kronig relation we have 
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e(g, uj) cu"^ — cu"^ + id 



(9) 



Within the PPA the self-energy can now be separated into a frequency independent exchange 
term and a correlation term 



(10) 



where 



Sex(fc, u;) = i I J^^Vc{q)Go{k + q,uj + uj), 
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and 
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GQ{k + q,uj + uj'). 



(12) 



The exchange energy Sex(fc,i^) as well as the correlation energy S^or within the full RPA 
theory has been calculated earlier by Hu and Das Sarmai. Using Eq. (^ in Eq. (|12D and 
performing a frequency integration, the correlation part becomes 



dq VMW p 
27r 2uJn 
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The real and the imaginary parts of the Scor are given by 



and 



dq 

ReT,coT{k,u) = P / ^9{(l) 
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\k + q\-kF)5{u;-ujq-^k+q)], (15) 



where g{q) = Vc{q)ujQ/ {2u!q) and P/ indicates the principle value integral. From the restric- 
tions on the integration region arising from various 6 and 6 functions, we see that ImScor is 
nonzero only for 



^q{kp — k) > uj and —Uq{kp + k) < u < if k < kp, 

^fc > and —ujqikp + k) < oj < —ojqik — kp) ii k > kp. 

Carrying out the integral over g, one obtains the imaginary part of the Ecor 
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ImScor(/^,tc;) = 
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\dn+{q+,)/dq\ \dn_{q_,)/dq\ 

where ^±{q) = ± cjg — ^^+5 and q±^i are zeros of fi±(g). From Eqs. (|T3|) and (|IB]) we 
know that Im[S(/c, cu)] as a function of uo has finite discontinuities aX uo = ±ujq{k + kp), 
whose magnitude can be calculated from Eq. (|l^. For example, we have the magnitude 
g{kp) /[dujq{kp) /dq ± kp/m] at u; = ±ujq{kp) for A; = and {l/2)g{2kp) /[dujg{2kp) /dq ± 
kp/m] at UJ = ±ujq{2kp) for k = kp. (See the numerically calculated values in Fig. 3(a) and 
(b).) A finite discontinuity in Im[S] gives rise to a logarithmic singularity in Re[S], which 
can be verified using the Kramers- Kronig relation. (See Fig. 3(a) and (b).) 

In order to determine quasiparticle excitation energies one must solve the Dyson 
equation^ which is given by 

uj + ij = ^{k) + i:{k,uj), (18) 

where /i is the chemical potential of the interacting electron gas, which is determined by 
setting k = kp and = in the above equation. Once the self-energy T,{k,uj) is known the 
single particle spectral function A{k,Lj) is readily calculated. A{k,u!) contains important 
dynamical information about the system and is given by 

2\linS(k,Lj)\ , 

A(k,uj) = ^ ^-^^ p. (19) 

[uj -^{k) -ReJ:{k,uj)f + [lmS{k,uj)f 

It satisfies the sum rule 

r^A{k,u) = l, (20) 



which we verify to be satisfied within less than a percent in our numerical calculations. 

Fig. 3 shows the calculated self-energies and spectral functions as a function of frequency 
u for = (band edge) and k = kp (Fermi energy). The complete RPA results^ (thin lines) 
are also shown for comparison with our PPA results. From the figures we can see that our 
PPA results are almost identical to the full RPA resultfi In both calculation the parameters 
corresponding to GaAs are used: m = O.OTmg (me is the free electron mass), cq = 12.9, 
€00 = 10.9, and culo = 36.8 meV. The well width of a = lOOA and the ID electron density of 
n = 0.56 X lO^cm"^, which corresponds to a Fermi energy Ep ~ 4.4meV and a dimensionless 
density parameter Vs = Ame'^ /Tckpeo = 1.4 with kp = 7m/2, are used in both calculation. In 
Figs. 3(a) and (b) the straight lines are given hj u — ^{k) — /i, and their intersections with 
Re[E] indicate the solutions to Dyson's equation and correspond to quasiparticle peaks. In 
the spectral function for A; = 0, we find two undamped quasiparticle peaks. The strength 
(27r X 0.37) of the regular quasiparticle (the first peak near u = 0) within PPA is slightly 
higher than the corresponding RPA result {2it x 0.33). The strength (27r x 0.31) of the 
second peak, the so-called plasmaron peak, is nearly the same in the PPA as that in the 
RPA (27r X 0.32). The low energy incoherent spectrum {Ep < uo < ujq{kp)) arising from 
the electron-hole continuum within RPA is transferred to the quasiparticle spectrum in the 
PPA, making the quasiparticle spectral weight slightly higher in the PPA than in the RPA. 
For k = kp we find that the quasiparticle-like peak at = is not a strict 5-function peak, 
which means that the system within PPA has no true long-lived quasiparticles. This result 
is qualitatively the same as the RPA result, implying that there can be no true quasiparticles 
in one dimension. As ^ 0, the dominant contribution to Im[S(/ci7, cu)] within RPA comes 
from the plasmon excitationil. Therefore, the behavior of the spectral function near u; = 
for k = kp shows exactly the same behavior for both the PPA and the RPA. 

In 2Di and 3Di§, the quantitative differences between the results of plasmon-pole ap- 
proximation and the random-phase approximation are comparable, and are considerably 
larger than what we find in our ID calculations. In our ID calculation, the agreement 
between the RPA and PPA self-energies is almost perfect. Since the ID electron-hole con- 

9 



tinuum is strongly suppressed by the severe phase restriction due to energy-momentum 
conservation, the ID plasmon is the dominant excitation which contributes to the electron 
self-energy, with the contribution from single-particle excitations being essentially negligibly 
small. In this paper we provide an easy method for calculating the effects of correlation 
on the single particle self-energy in one dimension. It should be fairly straightforward to 
extend the PPA self-energy calculation to more complicated experimentally relevant situa- 
tions, such as finite temperatures and multisubband occupancies, with reasonable confidence 
of obtaining quantitatively accurate results. This is the main significance of our work. 



III. ENERGY RELAXATION IN A QUANTUM WIRE 

In this section, we apply the plasmon-pole approximation to a coupled electron-LO- 
phonon system in a quantum wire and study hot electron energy relaxation^ through LO- 
phonon emission. Although this topic is of great importance by itseli0, the present purpose 
is to use it as an example to show the simplicity and the reasonable quantitative accuracy 
of the plasmon-pole approximation in calculating quantum wire electronic properties. We, 
therefore, refrain from discussing in details the hot electron relaxation phenomena0. 

When excess energy is supplied to an electron gas, the electrons go out of equilibrium 
with the underlying lattice, with the electron gas attaining an effective electron temperature 
T higher than the embedding lattice temperature T^. Such a hot electron gas loses energy to 
its surrounding in order to return to equilibrium with the lattice. In polar semiconductors 
such as GaAs, the most efficient energy relaxation process, except at very low electron 
temperatures, is through LO-phonon emission. When reabsorption of the emitted LO- 
phonons is ignored, the hot electron energy loss rate at zero lattice temperature is given 
bylli (we take = throughout, our results should be vahd for low values of T^): 

—unT{u)\M,\Hu,x"'\q,uj)luiD^''\q,u), (21) 

where nT{uj) is the Bose distribution factor at electron temperature T, and |Mqp is the 
Frohlich coupling matrixB. The phonon propagator in Eq. ( ^1]) is 
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(22) 



The last term in the denominator is the phonon self-energy correction due to many-body 
electron-phonon coupling, which broadens the phonon spectral function. The phonon mode 
couples to the plasmon excitation (the so-called plasmon-phonon coupling) as well as to 
particle-hole excitations, so that the renormalized phonon spectrum may be characterized 
as containing hybridized phonon-like and plasmon-like modes, and quasiparticle-like modes. 
The phonon-like mode has large spectral weight and high energy (~ ujlo), while the other 
modes have small spectral weights but have arbitrarily low energies. At high electron temper- 
atures (fcfiT ~ ojlo)i energy relaxation through emission of the phonon-like mode dominates 
because of its large spectral weight, while at low temperatures (fc^T ^ cj^o), energy re- 
laxation though emission of the plasmon- and quasiparticle-like modes dominates because 
of their low energies. The existence of the low energy modes enhances the energy loss rate 
at low temperatures since emission of bare phonon mode with frequency uolq is effectively 
frozen out when fc^T ^ ojlq- Our present objective is to compare the energy loss rates 
among these three cases: no many-body phonon-electron coupling, involving only the bare 
phonon mode; phonon-electron coupling in the plasmon-pole approximation, involving only 
the hybridized plasmon- and phonon-like modes, but no quasiparticle-like modes; and the 
phonon-electron coupling in the full RPA, involving all the modes. One can see that these 
are three increasingly sophisticated approximations to the phonon self-energy correction in 
Eq. (p2|) with the phonon self-energy correction completely neglected in the bare phonon 
case. 

With the phonon self-energy ignored, ImD(g, oj) becomes a single ^-function at = ojlo- 
Eq. (pH) then gives the energy loss rate as 



The characteristic of the bare phonon result is an approximate exponential temperature 
dependence Pq oc exp(— culo/^b^); which comes from the Bose factor nxiuji^o)- With the 




(23) 
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plasmon-pole approximation xppi Ini-D becomes a pair of 5-functions at = uo±, the fre- 
quencies of the hybridized plasmon-phonon modesEl. The energy relaxation rate is then 
given as 

Ppp = P+ + P_, (24) 

with 

P± = Y.^±M^± f''i'^l ~ "fJ \M,\\-2)\mx{q,^±). (25) 

where P± refer respectively to energy loss via upper (lower) hybrid plasmon-phonon modes. 
It should be noticed that the above expression is formally as simple as the corresponding 
bare phonon result given in Eq. (^), both involving a wavevector integral. 

The energy loss rates with no phonon renormalization, with phonon renormalization in 
the plasmon-pole approximation, and with phonon renormalization in the full RPA, are 
shown in Fig. ^. Two things need to be emphasized. The first is that the phonon renormal- 
ization enhances the energy loss rate by orders of magnitude at low temperatures, although 
its effect is negligible at high temperatures. The second is that the plasmon-pole approxi- 
mation gives an excellent description of the energy loss process, in the sense that its result 
agrees very well with the full RPA result. This example shows again that the plasmon-pole 
approximation can work remarkably well in a quantum wire system because of the increased 
dominance of plasmon excitation in one dimension. 

IV. SUMMARY 

The plasmon-pole approximation has been widely employed in three- and two- 
dimensional many-body electron systems. In this work, we discuss two specific applications 
of the plasmon-pole approximation to ID electrons in a quantum wire structure. Our re- 
sults suggest that the plasmon-pole approximation can work exceptionally well in calculating 
electronic many-body properties in a semiconductor quantum wire structure because of the 
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severe phase space restriction on single-particle electron-hole excitations in one-dimensional 
systems. We apply the plasmon-pole approximation to calculations of electron self-energy 
due to Coulomb interactions and hot electron energy relaxation rate via LO-phonon emis- 
sion, and find that our calculated PPA results agree extremely well with the results of the full 
RPA calculations. The agreement of the PPA results with the full RPA results is substan- 
tially better (in fact, essentially exact) in ID than in the corresponding 2D and 3D systems. 
Our results should influence future electronic calculations in semiconductor quantum wires 
where calculations may now safely ignore the full complications of the RPA and adapt the 
simple, intuitively appealing and quantitatively accurate PPA. 
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FIGURES 

FIG. 1. Excitation spectrum of a one-dimensional electron gas in the RPA. Particle- hole ex- 
citations arc confined within the phase space surrounded by the dottcd-linc ABODE. Plasmon 
excitation is represented by the solid-line. 

FIG. 2. Calculated RPA oscillator strengths of plasmon excitations of ID quantum wire and 2D 
quantum well electron systems. The input parameters are taken from GaAs-based materials: den- 
sity n = lO^cm"^ and lateral width a = b = 200A for the quantum wire; density n = 1.6x lO^-'^cm"^ 
for a zero thickness purely 2D quantum well. 

FIG. 3. (a), (b) Self-energy T,{k,Lo) and (c), (d) spectral function A(k,Lo) as functions of the 
frequency u for two fixed wave vectors A; = ((a) and (c)) and kp ((b) and (d)). Thick (thin) 
lines correspond to the PPA (RPA) results. The vertical lines in (c) represent (5-functions with the 
spectral weight given above the peaks. The straight lines in (a) and (b) are given by a; — ^(fc) — p,, 
and their intersections with Re[S] indicate the solutions to Dyson's equation and correspond to 
quasiparticle peaks. |ImE| is plotted instead of ImS for visual clarity. 

FIG. 4. Energy relaxation rates per electron as functions of electron temperature T. The dot- 
ted-line, dot-dashed-line, and solid-line are respectively the results with no phonon renormalization, 
with phonon renormalization in the plasmon-pole approximation, and phonon renormalization in 
the full RPA. The quantum wire has a carrier density n = 10^ cm""^ and wire widths a = b = 200A. 
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Figl Das Sarmo, Hwang, Zheng 
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